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Ab initio phonon calculations for the nonpolar reference 
structures of the (001), (110), and (111) PbZr0 3 /PbTi0 3 
[1:1] superlattices are presented. The unstable polar modes 
in the tetragonal (001) and (110) structures are confined in 
either the Ti- or the Zr-centered layers and display two-mode 
behavior, while in the cubic (111) case one-mode behavior 
is observed. Instabilities with pure oxygen character are ob- 
served in all three structures. The implications for the ferro- 
electric behavior and related properties are discussed. 



I. INTRODUCTION 

The soft-mode theory of ferroelectricsEI has established 
the relationship of the unstable polar modes of a high- 
symmetry nonpolar reference structure to the physics 
of both the ferroelectric transition and the large dielec- 
tric and piezoelectric coefficients. First-principles den- 
sity functional perturbation theory (DFPT),BB permits 
the direct quantitative investigation of such unstable 
modes and the accurate computation of the interatomic 
force constants, phonon frequencies and dynamical ma- 
trix eigenvectors of the cubic reference structures of fer- 
roelectric perovskite oxides such as BaTi03, PbTiC>3, 
SrTiC>3, and KNbC>3, has lead to a deeper understanding 
of the individual compounds and materials trends3l3. 

For technological applications requiring high- 
performance dielectrics, piezoelectrics, and ferroelectrics, 
the most promising candidates are not pure perovskite 
compounds, but solid solutions, Siich as Pb(Zr a; TL_ T2 ;)03 
(PZT)H, (Ba x Sr!_ x )Ti03 (BST)I and PMN-PTE3, and 
atomic-scale supedattices including BaTiOs/SrTiOs and 
KNb03/KTa03.llil In these systems, the unstable polar 
modes of their high-symmetry nonpolar reference struc- 
ture are similarly expected to be central to a predictive 
understanding of their properties. However, their chem- 
ical and structural complexity, requiring the use of large 
supercells, leads to a tremendous increase in the compu- 
tational demands of first-principles calculations. 

Accurate first-principles calculations can be performed 
only for short period superlattices and the simplest or- 
dered configurations of solid solutions. PZT is a natural 
choice for first-principles investigation. Not only is it al- 
ready in widespread use in piezoelectric transducer appli- 
cations, but its optimal composition, at the morphotropic 
phase boundary, is almost exactly x=0.5, which corre- 



sponds to ordered structures with supercells as small as 
double the primitive perovskite unit cell. Indeed, work 
on the (001) and (111) doubled supercells, dominates 
available first-principles results on PZT.EjTj While the 
atomic arrangements in real bulk PZT are of course much 
more disordered, the study of the lattice dynamics of 
ordered structures with [1:1] compositional modulation 
along various crystallographic directions provides (i) an 
accurate characterization of the phonons of ultrathin su- 
perlattices, (ii) a guide to understanding the vibrational 
properties of the disordered alloys, and (iii) a quantita- 
tive benchmark to develop and test ab initio based ap- 
proximations that can accurately describe the lattice dy- 
namics of disordered solid solutions and large scale het- 
erostructures, and are computationally more amenable 
than direct first-principles calculations. 

In this paper, we present first-principles calculations 
for the lattice dynamics of three ordered configurations of 
Pb(Zr .5Tio.5)0 3 : the [1:1] ultrathin superlattices (001), 
(110) and (111). In Section II, we give the details of 
the computations. In Section III, we present the relaxed 
nonpolar structures, phonon dispersions along selected 
directions, and detailed analysis of the unstable modes. 
Several interesting features not present in the pure end- 
point compounds, and the implications for the ferroelec- 
tric behavior and related properties, are identified and 
discussed in Section IV. 



II. COMPUTATIONAL METHOD 

We performed ab initio calculations using density func- 
tional theory (DFT) with the plane wave pseudopoten- 
tial method. The vibrational properties were computed 
using density functional perturbation theory (DFPT)O 
generalized to ultrasoft pseudopotentials.Ea All calcula- 
tions J&fre performed using the PWSCF and PHONON 
codeslij |— - 

We used Vanderbilt ultrasoft pseudopotentialsE3, 
treating as valence states the 3s, 3p, 3d, and 4s elec- 
trons of Ti, the 4s, Ap, Ad, and 5s electrons of Zr, the 5d, 
6s, and 6p electrons of Pb, and the 2s and 2p electrons of 
O. The inclusion of the semi-core electrons in the valence 
states is necessary for an accurate description of these ox- 
ides. The exchange and correlation energy is given within 
the local density approximation (LDA) using the param- 
eterization of Perdew and ZungerEl A kinetic energy 
cut-off of 35 Ry was used and the augmentation charges 
were expanded up to 350 Ry. For PZT(001), (110) and 
(111), the Brillouin zone (BZ) integration was performed 
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using 12 (663 tetragonal mesh), 9 (446 tetragonal mesh) 
and 6 (333 fee mesh) k-points in the irreducible BZ, re- 
spectively. 

For the prototypical cubic perovskite structure of the 
pure materials, PbZr0 3 (PZ) and PbTi0 3 (PT), we ob- 
tained the following values for the theoretical equilibrium 
lattice parameter: a(PZ)=7.77 a.u. and a(PT)=7.37 a.u. 
(experimental values 7.81 a.u. and 7.50 a.u.) and bulk 
modulus B(PZ)=170 GPa and B(PT)=202 GPa. These 
compare very well with results from previous LDA cal- 



culations, such as Ref. 22 



To compute the LO-TO splitting at T arising from 
the long-range Coulomb interactions present in the per- 
ovskite oxides, it is necessary to know the electronic 
dielectric tensor eoo and the Born effective charge ten- 
sors Z*, where s runs over the atoms in the unit cell. 
These can be computed within the framework of DFPT, 
although the current version of PHONON allows these 
computations only with norm-conserving pseudopoten- 
tials. For the calculation of and Z* we hajffi. therefore 
used Troullier and Martins pseudopotentialso with the 
same valence states as for the ultrasoft pseudopotentials. 
To achieve convergence with the norm-conserving pseu- 
dopotentials an energy cutoff of 80 Ry was needed. All 
the other parameters, such as lattice parameters, atomic 
positions, and k-point sampling were chosen to be the 
same as those optimized with the ultrasoft pseudopoten- 
tials. 



III. RESULTS 

The primary focus of our calculations is the eigenfre- 
quencies and eigenvectors of the unstable polar modes at 
r. The unstable eigenvectors of the analytic part of the 
dynamical matrix are of particular interest, since they 
generate symmetry-breaking energy-lowering distortions 
of the relaxed structures under conditions of zero macro- 
scopic electric field. Some direct connections between 
features of the superlattice structure and the unstable 
modes emerge. In general, we will see that if pure Ti-0 
chains are present in the ordered structure, the lowest 
frequency unstable modes are confined in the PT layer 
and are characterized by displacements parallel to these 
chains. Another connection is between modes of the su- 
perlattice structures and those of strained bulk PT and 
PZ. This is especially easy to see in the ideal structure, 
where the relevant strain is an isotropic expansion for PT 
and compression for PZ. For the noncubic (001) and (110) 
structures, the LO-TO splitting gives rise to a nontriv- 
ial anisotropy with changing direction of the vanishing q 
vector. Surprisingly, unstable LO modes are also found 
to occur. Throughout this section, symmetry analysis 
will be used to organize and present the results. 



A. Relaxed structures 

To find an appropriate high-symmetry reference struc- 
ture for PZT, we first placed the atoms as in the 
ideal cubic perovskite structure with lattice constant 
a(PZT)=7.57 a.u., obtained by averaging the theoreti- 
cal values in the pure endpoint compounds PZ and PT. 
This value is very close to the LDA value o£-Z.55 a.u. cal- 
culated for PZT(lll) by Fornari and Singh.E3 The Zr and 
Ti atoms are ordered on the B sites to establish the (001), 
(110) and (111) [1:1] superlattices. These superlattice 
structures will be referred to below as "ideal" PZT(001), 
PZT(llO) and PZT(lll). The "relaxed" superlattice 
structures are then obtained by allowing symmetry pre- 
serving relaxations of the internal atomic coordinates, 
holding the lattice parameters fixed. The three ordered 
supercells each contain ten atoms, with the Wyckoff la- 
bels and corresponding atomic positions indicating the 
free internal structural parameters given in Table |. 

The unit cell of PZT(001) is tetragonal, obtained by 
doubling the cubic perovskite cell along the [001] direc- 
tion (space group D^). The lattice parameters are fixed 
at a=a(PZT) and c=2a. Alternating Ti0 2 and Zr0 2 
planes are separated by PbO planes along the [001] di- 
rection, the oxygens in the Wyckoff positions 2f, 2e, and 
2g being coplanar with the Ti, Zr and Pb atoms, re- 
spectively. Pure Ti-0 and Zr-0 chains run along the 
[100] and [010] directions, perpendicular to the direc- 
tion of compositional modulation, [001]. The symme- 
try preserving relaxation of the internal atomic coordi- 
nates consists of displacements along z of the atoms in 
the PbO planes (Pb and 0(2g) atoms). In the ideal per- 
ovskite structure, the z values for Pb and 0(2g) are ex- 
actly zpb = zo = 0.25. After relaxation we find that 
the PbO planes shift almost rigidly towards the Ti0 2 
plane with Azpb = —0.20 a.u. and Azo = —0.19 a.u., 
resulting in a lowering of the energy by 291 meV per 
unit cell. The relaxed Ti-0(2g) (Zr-0(2g)) distance is 
smaller (larger) than the corresponding distance in the 
pure compound PT (PZ). Since the BOg octahedron can 
be regarded as the building block of the perovskite struc- 
ture, it is useful to compare the distortion of the TiOe 
and ZrOe octahedra in the supercell relative to the pure 
endpoint compounds. In the relaxed (001) structure the 
TiOe (ZrOe) octahedron is stretched (compressed) in the 
(001) plane, due to the larger (smaller) lattice parameter 
(a(PT)< a(PZT)< a(PZ)), and shrinks (expands) along 
the [001] direction of compositional modulation. This re- 
laxation relieves some of the superlattice epitaxial strain 
by restoring the local volumes towards their equilibrium 
values in the bulk endpoint compounds, though the local 
c/a ratio is significantly changed. 

The PZT(llO) structure is also tetragonal (space 
group T)\ h ) with lattice parameters a=\/2a(PZT) and 
c=a(PZT). Along the [001] direction, TiZr0 4 planes al- 
ternate with Pb 2 02- Within the TiZr04 planes the Ti 
and Zr atoms are placed in a 2D checkerboard arrange- 
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ment, each surrounded by four coplanar oxygen atoms in 
the 4j position. Pure Ti-0 and Zr-0 chains run along 
the [001] direction, perpendicular to the two equivalent 
directions of compositional modulation, [110] and [110]. 
The symmetry preserving relaxation is specified by a sin- 
gle parameter Ax, the displacement of the four 0(4j) 
atoms along the Ti-Zr lines in the TiZr04 planes. In 
the ideal perovskite structure the O atoms are exactly 
halfway between Ti and Zr, with x=0.25. After relax- 
ation we find Ax = —0.125 a.u., so that the four co- 
planar 0(4j) atoms surrounding the Ti (Zr) are 0.025 
a.u. closer (further) than in PT (PZ), resulting in a low- 
ering of the energy by 338 meV per unit cell. Compared 
to the pure endpoint compounds the TiOg (ZrOg) octa- 
hedron is stretched (compressed) along the [001] direc- 
tion, and shrinks (expands) in the (001) plane of com- 
positional modulation. As in PZT(001), this relaxation 
relieves some of the supcrlatticc epitaxial strain by restor- 
ing the local volumes towards their equilibrium values in 
the bulk endpoint compounds, though the local c/a ratio 
is significantly changed. 

The PZT(lll) structure is obtained by placing Ti and 
Zr atoms on the simple cubic lattice of B sites in the 3D 
checkerboard arrangement corresponding to the rocksalt 
structure (space group O^). The unit cell is face cen- 
tered cubic with lattice parameter a=2a(PZT). There 
are no pure Ti-0 or Zr-0 chains in this structure; in 
all B-0 chains the Ti and Zr alternate. The symmetry- 
preserving relaxation is specified by a single parameter, 
Ax, corresponding to a breathing of the oxygen octa- 
hedra around the Ti and Zr atoms. A negative value of 
Ax corresponds to a breathing-in (breathing-out) motion 
around the Ti (Zr) atoms. In the ideal perovskite struc- 
ture, the O atoms are exactly halfway between Ti and Zr, 
with x = 0.25. After relaxation we find that the oxygen 
atoms have moved closer to the Ti, with Ax — —0.115 
a.u., resulting in a lowering of the energy by 465 meV 
per unit cell. The superlattice TiOg and ZrOe oxygen 
octahedra are remarkably similar to those in the pure 
endpoint compounds. Thus, at least with respect to the 
BOe octahedra (the Pb atoms are fixed), the superlattice 
strain is completely relaxed. 



B. Z* and eoo 

The computed electronic dielectric tensor and the 
Born effective charge tensors Z* for the relaxed PZT(OOl) 
and PZT(llO) structures are given in Tables ^ and D. 
All matrices are diagonal except for ZL^.j in PZT (110) 
(note that in PZT(llO) the xy axes are rotated by 45 
degrees around z with respect to the original cubic axes), 
and therefore, except in that case, only the diagonal ele- 
ments are given. 

The anisotropy associated with the atomic ordering 
is generally quite small. Also, while Z* Zr and Z^ are 
quite different in the endpoint compounds, even when 



computed at the same lattice constant a(PZT), the ef- 
fective charges of Zr and Ti in the superlattices are very 
similar Jij As a result, the dielectric tensor and effective 
charges obtained by averaging the tensors for pure PT 
and pure PZ, also given in Tables |l| and [l| are very 
close to the corresponding quantities for the superlat- 
tices. In fact, the change in LO phonon frequencies ob- 
tained by using the average values of and Z* instead 
of the superlattice values is in each case less than 2%, 
comparable to the overall accuracy expected for the first- 
principles calculation. The average values might there- 
fore be expected to be a good approximation in other 
x=0.5 structures as well. The results that follow were 
computed using the average values. 



C. Phonon dispersion 

For the three relaxed (001), (110), and (111) PZT 
superlattices, the phonon dispersion along the direction 
of compositional modulation and the anisotropy at q = 
are shown in Figure 1. In this section, we give an 
overview of the results, with detailed discussion and anal- 
ysis to follow in Sections HID, E and F. 

The behavior of the phonon dispersion along the direc- 
tion of compositional modulation indicates the degree to 
which the phonon vibrations are confined to layers of a 
particular composition. This dispersion is quite flat for 
most of the modes in the (110) and, to an even greater de- 
gree, in the (001) superlattices. This means that the cor- 
responding normal modes are strongly localized within 
atomic planes perpendicular to the direction of composi- 
tional modulation. In contrast, most of the modes in the 
(111) superlattice display significant dispersion, indicat- 
ing that the normal mode amplitudes are spread through 
the supercell. This will be confirmed by the analysis of 
individual eigenmodes below. 

Due to the long-range character of the Coulomb force, 
the dynamical matrix displays non-analytic behaviour in 
the limit q^0, resulting in a dependence of the frequen- 
cies and eigenvectors on the direction of the vanishing 
wave vector, q. More precisely, in the long wavelength 
limit the dynamical matrix can be written as the sum of 
an analytic contribution, D T (q=0), and a non-analytic 
contribution that depends on the Born effective charge 
tensors Z* and the electronic dielectric tensor e„o: 



D ss ,( q ->o) = d;. 
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where D T (q=0) is the dynamical matrix computed with 
the boundary condition of zero macroscopic electric field 
(E=0), M s is the atomic mass of the atom s, and £1 is 
the volume of the unit cell.E-3i23 In the lefthand panels of 
the three plots in Figure 1, we show how the frequencies 
depend on the angle 9 between the direction of the van- 
ishing wave vector, q, and the direction of compositional 
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modulation. For a cubic structure, such as the (111) su- 
perlattice, the limiting frequencies do not depend upon 
the angle and can simply be identified as LO (longi- 
tudinal optic) and TO (transverse optic), though as in 
the perovskite structure, with multiple optic modes of 
the same symmetry, no one-to-one correspondence can 
be made between individual TO and LO modcsc3. For 
the tetragonal (001) and (110) superlattices, the limit- 
ing frequencies do depend on the angle 9. This gives rise 
to bands of LO frequencies which can be identified by 
their dispersion in Figure 1 and which will be discussed 
in detail in Section IIIF. 

All three superlattices display lattice instabilities; in 
fact, their phonon spectra contain entire branches of 
imaginary frequencies. The imaginary frequencies are in- 
dicated by the negative values in Figure 1 and the corre- 
sponding region of unstable modes is shaded in gray. The 
most unstable mode is observed at T in all three cases, 
consistent with our expectation that as for bulk PZT, 
the ground state of these superlattices should be ferro- 
electric. However, (111) is anomalous in that this dom- 
inant r mode is actually a nonpolar pure oxygen mode 
that will not produce a ferroelectric ground state, as will 
be discussed in more detail below. In the following, we 
focus on the behavior at the zone center, with particular 
attention to the unstable modes. 



D. Unstable eigenmodes of D T (q=0) 

Diagonalization of the analytic contribution to the dy- 
namical matrix, D T (q=0), gives the energies of the nor- 
mal modes in zero macroscopic field. These modes can 
be classified using the full symmetry of the point group 
of the space group, as will be described in detail be- 
low. The transverse IR-active phonons (TO) and the 
non-polar phonons are eigenmodes of D r (q=0) as well 
as of the full dynamical matrix, while the LO modes are 
mixed by the nonanalytic contribution, as will be seen in 
Section F. 

Since D T (q=0) has the full symmetry of the point 
group of the space group, we can classify its eigenmodes 
according to the irreducible representations (irreps) of 
D 4/l for PZT(OOl) and PZT(llO) and O^ for PZT(lll), 
following the nomenclature of Ref. 27. These labels spec- 
ify the polar character of the modes, which atoms are 
involved in the vibration, and the directions of the cor- 
responding atomic displacements. 

For the (001) supcrlattice, the eigenmodes are labelled 
by five irreps: 



6T y © 2Ti © 2P 
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where the coefficient gives the multiplicity of the irrep 
and the superscript indicates its dimension. The two- 
dimensional irreps F51 and L5 correspond to vibration in 
the xy-plane perpendicular to the direction of composi- 
tional modulation (_L). The others are one-dimensional 



irreps and correspond to vibrations along the direction 
of compositional modulation (||) . There are twelve IR- 
active polar TO modes (seven rg/(±) and five Li'(||) 
modes) in which all of the atoms participate. The re- 
maining T51 mode and Ty mode are the zero-frequency 
acoustic modes. All the other modes are non-polar: the 
two P3' modes involve only 0(2e) and 0(2f) z- vibrations, 
the two Ti modes involve only 0(2g) and Pb z- vibrations, 
and the two T$ modes involve only 0(2g) and Pb xy- 
vibrations. 

For the (110) superlattice, the eigenmodes are labelled 
by nine irreps: 



sr.; 



(2) m r( 2 ) 
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These irreps correspond to vibrations in the xy-plane of 
compositional modulation (||) or to vibrations in the per- 
pendicular z direction (_L). As in the (001) supercell there 
are twelve IR-active polar TO modes (seven (|| ) and 
five ri/(_L)) in which all of the atoms participate. The 
remaining mode and Fy mode are the zero-frequency 
acoustic modes. All the other modes are non-polar: the 
Ty mode involve only Pb z-vibrations, the L 4 ' and r§ 
modes involve only 0(4j) z-vibrations, and the Tx, L2, 
T3, r 4 modes also involve only 0(4j) vibrations but in 
the xy plane. 

For the (111) supercell the eigenmodes are labelled by 
nine irreps: 



5L 



(3) 



2r 



(3) 
25' 



r 



(3) 
25 



,(3) 
15' 
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There are four IR-active polar TO modes, Tis, which 
involve atomic displacements of all the atoms along x, 
y or z. The remaining L15 mode is the acoustic zero 
frequency mode. All the other modes are non-polar: the 
two T25' modes are the odd linear combination of the 
two Pb(x,y,z) and one oxygen displacement pattern, and 
the remaining T25, Tis', L12, and T\ are all unique pure 
oxygen modes. 

For both the ideal and relaxed structures of the three 
superlattices, the eigenfrequencies and symmetry labels 
of the unstable eigenmodes of D T (q=0) are given in Ta- 
ble IV . A comparison of the instabilities of the three ideal 



superlattice structures provides information on their sen- 
sitivity to the atomic configuration. In the ideal struc- 
tures the atomic positions are identical and the only dif- 
ference is in the arrangement of Zr and Ti atoms on the 
B sites. It can immediately be seen that the unstable 
eigenfrequencies are quite sensitive to the arrangement 
of atoms on the B site. In part, this reflects the differ- 
ent displacement patterns imposed by symmetry. Even 
if the Zr and Ti atoms were completely equivalent, these 
lists would differ because of the different q-points of the 
primitive Brillouin zone that are folded into the smaller 
cubic Brillouin zone of each supercell (L and X for (001), 
T and R for (111), and T and M for (110)), so only half 
of the modes (those from primitive T) would match for 
each pair. With the actual atoms in place the mixing 
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of the different "folded" mode combinations gives quite 
different results in the three superlattices. 

For each superlattice structure, the effect of atomic re- 
laxation has been investigated by comparing frequencies 
of the ideal structure with those of the relaxed struc- 
ture. In PZT(OOf) the shift in frequencies is substan- 
tial, and in some cases even leads to a reordering of the 
modes (for example, the unstable IV and T^> modes). 
The relaxations in the other two superlattice structures 
cause smaller changes. The larger shift in PZT(OOl) may 
be due to the larger atomic relaxations observed in this 
structure, which involve Pb displacements in addition to 
the oxygen displacements. 

The unstable polar TO modes will be described in de- 
tail in the next section. In addition, each of the three 
superlattice structures has one pure oxygen instability. 
The oxygen instability in PZT(OOl), the Ty mode, is 
mainly confined in the PZ layers and corresponds to the 
r 2 5 mode of bulk PZ. In the superlattice the PZ layers 
are under compressive strain (a(PZT)< a(PZ)) and the 
origin of this oxygen instability is in the pressure depen- 
dence of the mode of PZ which becomes unstable un- 
der compression. We note that this Ty oxygen instability 
is very sensitive to strain. In fact, the atomic relaxation 
alone causes a large shift of its eigenfrequency. The oxy- 
gen instability in PZT(llO), the r 2 mode, is a rotation of 
the oxygen octahedra around the [001] direction and cor- 
responds to an M mode of the primitive bulk cubic cell. 
In PZT(lll) the oxygen instability is the most unstable 
r 2 5 mode. It involves the rotation of the oxygen octahe- 
dra around the < 1 1 1 > directions and corresponds to the 
bulk R25 mode. This unstable mode contributes to the 
orthorhombic antifcrroelectric ground state of pure PZ 
and is responsible for the tetragonal antifcrrodistortive 
ground state of SrTiC>3. Comparing the instabilities of 
PZT(lll) with those of PZT(OOl) and PZT(llO), we see 
that the frequency of the oxygen instability is similar in 
all three systems. However, in PZT(lll) the most un- 
stable polar mode has been significantly shifted upwards, 
which we attribute to the lack of intact Ti-0 chains in 
any direction. This leaves the oxygen instability as the 
dominant unstable mode. 



E. Confined TO polar instabilities 

Each of the three superlattice structures possesses TO 
polar instabilities, indicated in boldfaced type in Ta- 
ble IV, that can generate candidate ferroelectric struc- 
tures. PZT(OOl) and PZT(llO) exhibit "two-mode be- 



havior," that is, there are two types of such unstable 
modes, one with atomic displacements mainly confined in 
the PT layer and the other in the PZ layer. The atomic 
displacements of the two most unstable TO modes of re- 
laxed PZT (001), T 5 / and IV, are shown in Figure 2. 
The Ts> is a PT-like mode polarized parallel to the in- 
terfaces; it is confined in the PT layer and is remarkably 



similar to the unstable TO mode of bulk PT (also shown 
in Figure 2 for comparison). The Ty mode, on the other 
hand, is a PZ-like mode polarized perpendicular to the 
interfaces; it is mainly confined in the PZ layer and is 
remarkably similar to the unstable TO mode in bulk PZ. 
The confinement of these modes in the PT and PZ layers 
is consistent with their flat dispersion for q perpendicular 
to the interfaces. 

The atomic arrangements and projections in the (110) 
and (111) structures do not lend themselves well to repre- 
sentation in a figure analogous to Figure 2. A convenient 
way of characterizing the confinement of the polar insta- 
bilities in all three superlattice structures is to decompose 
the mode effective charge, Z, into two contributions: one 
from the Zr-centered cell and one from the Ti-centered 
cell: Z = Zz r + Zxi- These local contributions to the 
mode effective charge, which will be referred to as "local 
effective charges," are defined as 



^Ti(Zr) 



= > W,Z. 



E 

s ,0 



s^af3;s u P,s 
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where the sum over the atoms, s, runs only over the 
atoms in the Ti (Zr) centered cell and the weighting fac- 
tor w s is w=l for Ti (Zr), w—0.5 for the 6 oxygen atoms 
and t«=0.125 for the 8 Pb atoms. In the case that the 
modes are confined, we expect the local effective charge 
to have its largest value in the corresponding layer. 

The local effective charges for PZT(001) have been in- 
cluded in Figure 2. For the PT-like mode, r 5 >, Z Ti =7.29 
is more than three times larger than Zz r , indicating a 
high degree of confinement. Moreover, it is very close to 
the mode effective charge of the unstable TO mode in PT 
(at the equilibrium lattice constant), Z=7.57, confirming 
the local similarity to the bulk PT unstable mode. For 
the PZ-like mode, Ti> , the larger contribution is from the 
Zr-centered cell, Zz r =4.93, which is likewise very close to 
the mode effective charge of the unstable TO mode in PZ, 
Z=4.86. 

The local effective charges of the unstable TO modes of 
all three superlattice structures are given in Table [V| By 
symmetry, the induced polarizations are parallel to the 
direction of atomic displacement for each mode, also in- 
cluded in the table. For almost every unstable TO mode 
of the (001) and (110) superlattices, either Zti or Zz T 
strongly dominates, indicating that the polar instabili- 
ties are mainly confined either to the PT or to the PZ 
layers, respectively. The confined character of the modes 
is seen to be rather insensitive to the relaxation. Just 
as for the (001) case discussed above, the most unstable 
mode for PZT(llO) is a PT-like mode with polarization 
aligned in the PT plane (perpendicular to the direction 
of compositional modulation, _L) followed by a PZ-like 
mode with polarization parallel to the direction of com- 
positional modulation, ||. 

In contrast to the two-mode behavior of PZT(001) and 
PZT(llO), PZT(lll) exhibits one-mode behavior. For its 
single unstable TO mode, the local effective charges are 
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almost identical and the polar instability is distributed 
throughout the supercell. The different behavior of PZT 
(111) is a consequence of the fee symmetry; in fact, the 
displacements of the six oxygens in the supercell are 
symmetry-related, so that any mode in which oxygen dis- 
placements are allowed will in general involve the whole 
supercell, preventing the confinement to occur. 

There are two distinct factors that could produce the 
observed mode confinement: the difference between the 
atomic mass of the two species occupying the B sites (Ti 
and Zr) or the dependence of the interatomic force con- 
stants upon atomic configuration (i.e. the atomic force 
constants may depend upon the chemical identity of ei- 
ther of the two interacting atoms and upon the surround- 
ing atomic configuration). In some, semiconductor het- 
erostructures, such as AlAs/GaAsc3'Ea, the interatomic 
force constants (IFCs) are to a good approximation un- 
affected by the substitution, Al/As, and the observed 
character of the heterostructure modes is the result of 
the difference in mass. One might hope that the isoelec- 
tronic substitution Ti/Zr could be similarly understood. 
To check this, the eigenfrequencies and local effective 
charges of the dynamical matrices in which the masses of 
Ti and Zr have been replaced by their average are shown 
in Table VI. Changing the value of the masses results 



in a frequency shift that is larger for those modes that 
involve larger displacements of the Ti or Zr ions. In par- 
ticular the PT-like modes of (001) and (110) are shifted 
to higher frequency. The PZ-like modes are almost un- 
affected since the atomic displacements of Zr and Ti are 
almost zero in these modes. By the analysis of the local 
effective charges in Table VI it is evident that the mode 
confinement is still present when Ti and Zr have the same 
mass. It must therefore be that the IFCs depend strongly 
upon atomic configuration, with interactions involving Zr 
atoms significantly differing from their Ti counterparts. 
The important consequence of this observation is that 
approximations for the IFCs that neglect the effect of 
compositional modulation, with the most widely used be- 
ing the virtual crystal approximation, cannot adequately 
capture the character of the instabilities in PZT. 



F. LO-TO splitting and anisotropy at q=0 

The effect of adding the non-analytic term to the dy- 
namical matrix is to mix the IR-active modes that have a 
non- vanishing longitudinal component of the polarization 
thereby generating longitudinal IR-active modes (LO). 
This results in a splitting of the LO and TO frequencies 
for vanishing q vector that, in general, depends upon the 
direction (as noted in the angular dispersions of figure 

The angular dispersions of the LO modes have been 
already displayed in Figure [v|. For a more detailed dis- 
cussion, and to assist in the interpretation of the figure, 
the frequencies of all the TO and LO modes of the three 



relaxed supercells are given in Table VII with their sym- 
metry labels and mode effective charges. For the (001) 
and (110) supcrlattices the LO frequencies are direction 
dependent and we give their values for q parallel and per- 
pendicular to the direction of compositional modulation 
in each supercell, corresponding to 9 = and 9 = 90 deg 
in the angular dispersion shown in Figure 1, respectively. 
For the cubic (111) superlattice, symmetry considera- 
tions show that the LO frequencies do not depend upon 

q- 

Examination of the table shows that the (001) and 
(110) superlattices have unstable LO modes, a feature 
not present in the pure compounds. These modes are 
not confined, showing significant relative motion of the 
cations and anions in both the Zr and Ti centered subcells 
of the supercell. However, the induced polarizations in 
the two subcells are opposite in direction, and resulting 
near cancellation makes the unstable LO mode effective 
charges much smaller than those associated with the TO 
instabilities. 



IV. DISCUSSION 

As for the pure endpoint compounds, the unstable 
TO phonons are expected to be valuable indications of 
the existence and nature of a ferroelectric ground state 
for the PZT superlattices here considered, as they gen- 
erate energy-lowering distortions with nonzero sponta- 
neous polarization P s under conditions of zero macro- 
scopic electric field (E = 0). 

With dominant TO instabilities, PZT(001) and 
PZT(llO) are virtually guaranteed to have ferroelectric 
ground states. In previous first-principles studies, atten- 
tion has focused on ferroelectric structures of PZT(001) 
with polarization along the [001] direction. The two- 
mode behavior discussed here suggests that the true 
ground state could be both more complicated and more 
interesting. Freezing in the unstable TO phonons of the 
(001) superlattice give rise to a macroscopic polariza- 
tion that is parallel to the interfaces in the PT layers for 
IV and perpendicular to the interfaces in the PZ layers 
for Ty. Depending on the anharmonic coupling between 
these modes and their coupling to homogeneous strain, 
either or both of these modes could freeze in to give a 
ground state structure with polarization along the [001] 
direction (the c phase), along [100] (the a phase), along 
[110] (the aa phase), in the (010) plane (the ac phase), 
or in the (110) plane (the designations of the phases cor- 
respond to those used in Ref. ^ in a different context). 

The choice of equilbrium structure type from the list 
above should also be quite sensitive to the mechanical 
and electrical boundary conditions. Applied stresses, epi- 
taxial constraints and applied electric fields will result 
in relative shifts of the two modes and changes in their 
coupling, and could result in transitions from one phase 
to another. In particular, this system might provide a 
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microscopic realization of a large-strain electrostrictive 
actuator as recently suggested by BhattacharyaH 

In contrast, PZT(lll) has a dominant nonpolar insta- 
bility, with the ferroelectric TO mode slightly higher in 
frequency. The ground state has been analyzed in de- 
tail by Fornari and SinghEj. They observe, consistent 
with our results, that the ferroelectric and rotational dis- 
tortions are quite competitive in energy and that this 
competition is very sensitive to strain. 

The (001) and (110) supercells display LO as well as 
TO instabilities, a feature not present in the pure com- 
pounds. This presents the interesting possibility of a per- 
pendicularly polarized ferroelectric ground state struc- 
ture for a freestanding PZT film or nanoparticle. The 
energy of such a structure, with a depolarizing field 
E = — 47rP, would in general be prohibitively high, cor- 
responding to the stiffness of the LO mode. However, 
with an unstable LO mode, there would be a net energy 
gain despite the depolarizing field. Since Z* is rather low 
for these modes, the spontaneous polarization would pre- 
sumably be quite small, though still possibly observable 
or even useful. Whether this instability can be made to 
yield an actual ground state through application of ap- 
propriate mechanical boundary conditions remains to be 
investigated. 

Aside from the physical implications of the calculated 
phonons, these results can also be used as a test of ap- 
proximate schemes for obtaining the phonons of these 
and larger ordered superstructures. It is clear that while 
the virtual crystal approximation has been seen to work 
well for PZT(lll), the requisite one-mode behavior of 
this structure is actually quite atypical. Generalizing 
from the lower symmetry PZT(OOl) and (110) structures, 
we expect that larger superstructures, with even lower 
symmetry, will exhibit multi-mode behavior. A key in- 
gredient in reproducing this behavior is that approxima- 
tions should allow for force constants involving Ti dis- 
placements to be different from the analogous force con- 
stants involving Zr displacements. The development and 
application of a suitable method, which was proposed in 
Ref. ||, will be the subject of a future publication.^ 

V. SUMMARY 

We have presented an ab initio study of the lattice 
dynamics of the high-symmetry nonpolar reference struc- 
tures of three ordered configurations of PZT: the three ul- 
trathin PT/PZ [1:1] superlattices (001), (110) and (111). 
In particular we focused on the zone-center lattice in- 
stabilities that can generate symmetry-breaking energy- 
lowering distortions and eventually lead to a phase tran- 
sition. In PZT(001) and (110), these show a two-mode 
behavior that may have interesting consequences for fer- 
roelectricity in these systems. Some direct connections 
between features of the superlattice structure and the 
unstable modes emerge. In general, if pure Ti-0 chains 



are present in the ordered structure, the lowest frequency 
unstable modes are confined in the PT layer and are 
characterized by displacements parallel to these chains. 
Another connection is between modes of the superlattice 
structures and those of strained bulk PT and PZ. Finally, 
we have verified that the character of the modes is not 
simply due to the small difference between the atomic 
mass of the intermixing cations, Ti and Zr, as is the case 
in some semiconductors, such as (Al,Ga)As, but that the 
effects of compositional modulation on the IFCs play an 
important role. Consequently, simple approximations for 
the IFC that neglect the effects of alloying, such as the 
virtual crystal approximation cannot provide an accurate 
description of the lattice dynamics in these systems. 
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FIG. 2. Atomic displacements for the two most unstable 
TO polar modes in PZT(OOl), Fg and r' 1; and the TOl un- 
stable modes of bulk PbTi03 and PbZr03. The projection 
in the (100) plane is shown. The different atomic species are 
represented as follows: small black circles for Pb, large white 
circles for O, small white circles for Ti, and small white cir- 
cles with a cross for Zr. Note that the Zr displacement is 
close to zero in all modes shown and the atomic vibration in 
the Zr-centered cell is dominated by the displacements of the 
oxygen atoms. 



TABLE I. Wyckoff labels and atomic positions, in reduced 
coordinates, for the three superlattice structures. For the 
(111) superlattice structure, the reduced coordinates are in 
terms of the lattice vectors of the conventional fee unit cell. 
For Wyckoff positions with multiplicity greater than one, the 
first atom in the list will be referred to as the "representative 
atom." 



PZT fOOll 


Pb(2h) 


(0 5 5 +zpul 




Ti(la) 


(0 0) 




Zr(lb) 


fO 51 




0(2e) 


fO ±zn) 




0(2f) 


(0 0.5 0), (0.5 0) 




0(2e) 


(0 0.5 0.5), (0.5 0.5) 


PZT (110) 


Pb(2e) 


(0.5 0.5), (0.5 0.5) 




Ti(la) 


(0 0) 




Zr(lc) 


(0.5 0.5 0) 




0(4j) 


(±x ±x 0) 




O(lb) 


(0 0.5) 




O(ld) 


(0.5 0.5 0.5) 


PZT (111) 


Pb(2c) 


±(0.25 0.25 0.25) 




Ti(la) 


(0 0) 




Zr(lb) 


(0.5 0.5 0.5) 




0(6e) 


(±x 0), (0 ±x 0), (0 ±x) 



TABLE II. Electronic dielectric tensor and Born effec- 
tive charge tensors for representative atoms of the relaxed 
PZT(001) structure. Only the diagonal elements are given; 
the off-diagonal elements are zero by symmetry. The dielec- 
tric and Born effective charge tensors obtained by averaging 
the tensors computed for pure PT and pure PZ (at the av- 
erage lattice parameter a(PZT)) are given on the right hand 
side for comparison. 







PZT(001) 






average 




Pb 


3.91 


3.91 


3.87 


3.90 


3.90 


3.90 


Ti 


6.51 


6.51 


6.36 


6.47 


6.47 


6.47 


Zr 


6.33 


6.33 


6.40 


6.47 


6.47 


6.47 


0(2g) 


-2.54 


-2.54 


-5.20 


-2.54 


-2.54 


-5.29 


0(2f) 


-2.77 


-5.33 


-2.28 


-2.54 


-5.29 


-2.54 


0(2e) 


-2.36 


-5.12 


-2.78 


-2.54 


-5.29 


-2.54 


too 


7.31 


7.31 


7.30 


7.61 


7.61 


7.61 
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TABLE III. Electronic dielectric tensor and Born effective 
charge tensors for the relaxed PZT(llO) structure. Only the 
diagonal elements are given, except for the 0(4j) atoms; the 
off-diagonal elements are zero by symmetry. The dielectric 
tensor and Born effective charges obtained by averaging the 
tensors computed for pure PT and pure PZ (at the average 
lattice parameter a(PZT)) are given on the right hand side 
for comparison. Note that in the (110) supercell the axis are 
rotated by 45 degrees around z with respect to the original 
cubic axes. 



PZT(llO) 


Average 


Pb 


4.24 3.62 3.93 


3.90 3.90 3.90 


Ti 


6.47 6.47 6.80 


6.47 6.47 6.47 


Zr 


6.47 6.47 6.24 


6.47 6.47 6.47 


O(lb) 


-2.49 -2.49 -5.48 


-2.54 -2.54 -5.29 


O(ld) 


-2.61 -2.61 -5.07 


-2.54 -2.54 -5.29 




/-3.90 -1.36 \ 


/-3.92 -1.38 \ 


0(4j) 


-1.36 -3.90 


-1.38 -3.92 




V -2.58 / 


V -2.54 / 




7.31 7.31 7.41 


7.61 7.61 7.61 



TABLE IV. Frequencies and symmetry labels of the un- 
stable eigenmodes of D T (q=0) for the relaxed and ideal 
PZT(001), (110) and (111) structures. The values in bold- 
faced type correspond to IR-active TO modes, the others be- 
ing nonpolar modes. In PZT(001) and PZT(llO) the corre- 
sponding atomic displacements are either parallel (||) or per- 
pendicular (_L) to the direction of compositional modulation, 
as indicated. Frequencies are in cm -1 . 







(001) 
relaxed 


ideal 






(110) 
relaxed 


ideal 




(in) 

relaxed ideal 


r S ' 


_L 


-185 


-235 


IV 


_L 


-231 


-240 


F25 


-144 -155 


IV 


II 


-177 


-151 


r 5 < 


II 


-148 


-158 


ris 


-110 -115 


r 3 ' 


II 


-153 


-72 


r 2 


II 


-134 


-147 






r S ' 


_L 


-73 


-114 


r 5 < 


II 


-76 


-71 






r B 


± 


-17 


-22 


rv 


_L 


-37 


-67 







TABLE V. The mode effective charge tensors, Z, associ- 
ated with the unstable IR-active TO modes are decomposed 
into two contributions, one from the Zr-centered cell and one 
from the Ti-centered cell: Z = Z Zr + Z T i- Results are shown 
for the relaxed and ideal PZT(001), (110) and (111) struc- 
tures. The mode effective charges for the unstable TO modes 
of bulk PbZr03 and PbTiOa are given for comparison, both 
at the equilibrium volume and strained to the average volume 
of the 50% alloy, V = a 3 (PZT). 



superlattice 






relaxed 


ideal 
Zzi Zti 




r 5 / 




2.21 


7.29 


1.57 


8.11 




r v 


(II) 


4.93 


2.99 


4.44 


3.71 




r 5 < 


(-L) 


3.38 


-1.56 


3.76 


-0.65 


(110) 


IV 


(-L) 


2.26 


7.73 


2.19 


7.87 




r 5 < 


(II) 


4.38 


2.81 


4.11 


2.91 




r 5 < 


(II) 


0.13 


3.03 


0.60 


3.42 




r v 


(-L) 


2.74 


0.04 


3.10 


-0.06 


(111) 


r 15 




3.81 


3.81 


3.71 


4.00 


bulk 






equilibrium 


strained 


PbZr0 3 
PbTiOa 






4.86 


7.57 


4.87 


8.00 



TABLE VI. The frequencies and local effective charges 
obtained for the relaxed superlattice structures when the Ti 
and Zr atoms are assumed to have the same atomic mass: 
M T i=Mz r =(M T i+Mz r )/2. 



superlattice 






Zzr 


Zt\ 


w(cm 


(001) 


r B ' 


(-L) 


2.32 


6.76 


-177 




IV 


(II) 


4.98 


2.99 


-178 




IV 


(-L) 


3.36 


-1.68 


-70 


(110) 


IV 


(-L) 


2.42 


7.12 


-214 




r B < 


(II) 


4.39 


2.80 


-148 




IV 


(II) 


0.11 


3.02 


-76 




IV 


(-L) 


2.82 


0.00 


-35 


(111) 


Fl5 




3.80 


3.82 


-110 
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TABLE VII. Frequencies and symmetry labels of the TO and LO zone-center eigenmodes. For PZT(OOl) and PZT(llO) the 
LO frequencies are given for the two directions qx and q\\ that are, respectively, perpendicular and parallel to the direction of 
compositional modulation. For PZT(lll), the LO frequencies are independent of the direction of q. For all modes, the mode 
effective charges are given in parentheses. 





TO 


(001) 
LO [q±] 


LO [q,,] 




TO 


(110) 
LO [qx] 


LO [qy] 




(111) 
TO 


LO 


F r / 


-185 (9.50) 


-82 (0.32) 




IV 


-231 


(9.99) 


-76 (1-16) 




Tis 


-110 (7.62) 


95 (1.34) 




-177 (7.92) 




94 (1.12) 


r 6 ' 


-148 


(7.19) 




-85 (0.39) 


Tis 


159 (9.27) 


330 (3.32) 




-73 (1.82) 


86 (0.25) 






-76 


(3.16) 




95 (0.72) 


Tis 


354 (3.27) 


396 (4.23) 


Fr/ 


88 (0.59) 


108 (0.95) 




IV 


-37 


(2.78) 


128 (2.11) 






534 (5.65) 


693 (9.90) 


IV 


110 (3.53) 




225 (1.73) 


IV 


98 


(0.68) 




112 (1.00) 








r S ' 


127 (4.03) 


193 (0.22) 




IV 


161 


(4.00) 




194 (0.54) 










194 (0.62) 


373 (5.31) 




IV 


161 


(4.78) 


370 (5.49) 










IV 


243 (1.83) 




306 (4.72) 


r r / 


206 


(2.96) 




320 (3.83) 








r r / 


440 (1.73) 


466 (3.80) 




r r / 


378 


(4.18) 




415 (3.53) 








IV 


465 (1.48) 




467 (0.72) 


IV 


417 


(1.12) 


434 (3.54) 










IV 


545 (6.86) 




716 (10.11) 


IV 


532 


(5.69) 




693 (9.94) 








IV 


658 (5.23) 


757 (8.56) 




IV 


661 


(5.35) 


760 (8.49) 
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